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We have investigated the quantum admittance of an ultra- 
smaU tunnel junction with arbitrary tunneling strength under 
an electromagnetic environment. Using the functional inte- 
gral approach a close analytical expression of the quantum 
admittance is derived for a general electromagnetic environ- 
ment. We then consider a specific controllable environment 
where a resistance is connected in series with the tunneling 
junction, for which we derived the dc quantum conductance 
from the zero frequency limit of the imaginary part of the 
quantum admittance. For such electromagnetic environment 
the dc conductance has been investigated in recent experi- 
ments, and our numerical results agree quantitatively very 
well with the measurements. Our complete numerical results 
for the entire range of junction conductance and electromag- 
netic environmental conductance confirmed the few existing 
theoretical conclusions. 

FAGS numbers: 73.23.Hk, 85.30.Wx, 73.40. Rw 



I. INTRODUCTION 

In semiclassical theory of Coulomb blockadeBi, elec- 
tron tunneling is controlled by the charging energy which 
is determined by the junction capacitance. While the ex- 
istence of this elementary chaining effect is clearly seen 
in multijunction configurationsO, in a single junction, be- 
cause of the stray capacitance, only partial blockade has 
been observed.Q The question is then, including the elec- 
tromagnetic environment, what is the effective junction 
capacitance? The work of Devoret et alu has shown that 
external electrical circuit has strong impact on the ef- 
fect of Coulomb blockade in a single junction, and un- 
der standard experimental condition where the external 
resistance is much smaller than the resistance quantum 
RK=h/e'^, the Coulomb gap is totally smeared out at 
low voltages and at low temperature. A similar result 
was also obtained by Girvin et aln In both works, as 
well as in most existing studies on charging ultrasmall 
tunnel junctions, the tunnel Hamiltonian was treated 
perturbativelyuu using the Fermi golden rule. This low- 
est order perturbative approach is valid for weak tunnel- 
ing strength. 

With increasing tunneling strength, one must go be- 
yond the Fermi golden rule to include higher order tun- 
neling processes. The conventional perturbative calcula- 



tion is very tedious, and has been used to study tiiQ-ata- 
tistical properties of a single electron box (SEB).Bll3 In 
a SEB an island is formed between a tunnel junction and 
a gate capacitor. This is the simplest structure exhibit- 
ing the charging effect. Even for such simplest system, 
it is very difficult to derive the partition function beyond 
the second order of the dimensionless tunneLepnductance 
Rk/Rt, where Rt is the tunnel resistance.E3 Therefore, 
an entirely different approach is required to investigate 
the physics associated to strong tunneling strength. 

The functional integral method has rbycn used to study 
an open, undriven tunnel junctionO^Ej As an extension 
of this work, the relevant properties of a SEB have been 
investigated with tiierM,th integral represeatation of the 
partition function.llj II2 It has been provedlj that a for- 
mal perturbation expansion of the partition function re- 
sults in exactly the same form as that derived with the 
path integral method, provided that the channel num- 
ber of the tunnel junction is very large, which is the case 
for metallic tunnel junctions. Furthermore, based on the 
functional integral formulism, a Monte-Carlo simulation 
can be constructed for the entire regime from weak to 
strong tunneling strength. The Monte-Carlo simulation 
has been performed for the SEB and a very smooth inter- 
polation between the weak tunneling and the strong tun- 
neling regime has been obtain.113 Therefore, functional 
integral is a reliable powerful tool to study ultrasmall 
metallic tunnel junctions with arbitrary tunneling rate. 

Since most experimental investigations on ultrasmall 
tunnel junctions measure the current response to applied 
voltage, a relevant theoretical study must take into ac- 
count the influence of the electromagnetic environment. 
This is the goal of the present work which will use the 
function integral approach. In Sec. II we introduce the 
model Hamiltonian of a single tunnel junction under an 
electromagnetic environment in a general form, and per- 
form a theoretical analysis for arbitrary tunneling rate. 
In Sec. Ill we continue to derive the quantum admit- 
tance of such system in a close analytical form. Then, 
in Sec. IV we restrict ourselves to the simple controllable 
electromagnetic environment which allows us to conmaxe 
our theoretical results with the recent measurements Ed'EJ 
For this case we only need the dc conductance which 
is the zero frequency limit of the imaginary part of the 
quantum admittance. The numerical results will be.pr|P-. 
sented in Sec. V which agree well with experimcntsEB'EJ 
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Finally, in Sec. VI we give a concluding remark. 

II. THEORY FOR MODEL HAMILTONIAN 

The Hamiltonian HQ=Hs+Ht of an open, undriven 
tunnel junction consists of two parts. The first one 

Hs=H, + H, (1) 

represents the system in the absence of tunneling, where 

H, = (2) 

is the Coulomb charging energy when a charge Q is added 
to the junction with capacitance C. The electrodes are 
modelled as free quasi-particle systems described by 



(3) 



where the channel index a includes both the transverse 
and the spin quantum numbers. The longitudinal wave 
vectors in the left electrod is labelled by k, and in the 
right electrod by q. c^^ (or c^^) is the creation operator 
of an electron with energy eua (or Eqa) in the left (or 
right) electrode. The second part of H describes the 
tunneling processes of electrons 



Ht 



kqa 



-iif) H- H.c. 



(4) 



where t^qa is the transition amplitude of an electron tun- 
neling from the state qa in the right electrode to the state 
k(7 in the left one. During the tunneling the channel index 
(T is invariant. Here we have neglected the tunneling time, 
which is approprite for metallic tunnel junctions, |0|ther- 
wise the tunnel Hamiltonian must be modified. Ea^Eil The 
operator ip is conjugate to Q via the commutation rela- 
tion 

[(p,Q] = -ie. 

To avoid possible ambiguity in our theoretical analysis, 
before introducing the electromagnetice environment, let 
us first write down the functiouaLintegral representation 
of the partition function of iJotJ't3 



Dipe 



-SoMt)] 



In the action Sq [(p] = Sc [(p] + St [(p] , 

S'r[</'1 — / dr— — 
Jo 425, 



(5) 



(6) 



describes the kinetic part of the system, where Ec = 
e^/2C is the elementary charging energy, and 

St[ip] = - / dr dr'atir - r') cos[(^(r) - ^{t')] (7) 
Jo Jo 



is resulted from the influence of tunneling processes. 
With the commonly used approximation that the transi- 
tion amplitude tkqa in (|^) is a constant tkqa—t, the junc- 
tion resistance Rt can be calculated as 

Rt = {A7r^\t\^piPrN)-\ 

where pi (or p,.) is the density of state in the left (or 
right) electrode, and N is the number of channels. In 
this case the damping kernel atir) can be expressed as 



at 



1 



4/32 sin^TTT/P)' 



(8) 



where the coefficient at—Rx/RT is the dimensionless 
junction conductance. The Matsubara frequencies loi = 
2t:1 I P which will be relevant to our future analysis and 
calculation lies in the region much smaller than the band- 
width D of the metal. For such frequencies the Fourier 
transform of atir) reduces to the simple form 



at{u;i) 



47r 



(9) 



The electromagnetic environment of the single tun- 
nel junction represented by Hq^ caused by the on-chip 
leads and pads etc., will be modelled by a transmis- 
sion line,El w^fj^H be described by a Caldeira-Leggett 
Hamiltonian 



2Cn 2e2 



(10) 



for an infinite series of LC circuits. In iJcx the first part is 
the charging energy of all capacitors and the second part 
is the magnetic energy of the corresponding inductors. 
The eigenfrequency of the nth LC circuit is given by i^n — 
1 / VLnCn, and the spectral density 



5{v - l^n) 



(11) 



is determined by the external impedance. The operators 
ipn and qn obey the commutation relation 



{^m 7 Qn 



-ieSr, 



The electromagnetic environment H^x is coupled to the 
tunnel junction Hq bilinearly 



n— 1 



E 



(12) 



The second term at the right hand side of the above 
equation is a physical counter term similar to the one 
introduced in the formulation of the dissipative quantum 
Brownian motion.EII 

The partition function Z of the total Hamiltonian 
H—Ho+Hex+Hint can be similarly obtained in the func- 
tional integral representation as 
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n=l 

exp \ ~So[tp] ' 



(13) 



dT 



2e2 



The path integrals over the bath modes are Gaussian, 
and thus can be evaluated exactly. From the condition 
of vanishing first variation of the action with respect to 
ifn, we obtain the classical equation of motion for the 
phase 

The solution to this linear differential equation with the 
boundary conditions 

fn{f3) = <y3„(0) = (fnO 

is 

(ci)/ N sinhK(/3-T)] 



93riO H / dr' smh{i^nT')(p{T') 

Jo 

sinh(i^„T) 
sinh(t/„/3) 



1 



dT'sinh[(/3-i.„)TXr') 



(14) 



The corresponding action is then calculated as 

= (15) 

(cV) 

where Siio depends only on the boundary value (p„Oi and 

S'exM = dr dr'ae^ir - t')[(p{t) - ip{T')f. 

^ Jo Jo 

(16) 

The damping kernel of the electromagnetic environment 
ctcx{T) can be expressed in a Fourier series 



l — ~OG 

with the Fourier coefficients 



acx(wi) 



dv J{v) Lof 

TT 1/ Ujf + 



(17) 



(18) 



related to the spectral density J(v). oicyA^oi) is thus pro- 
portional to the Fourier tramsform of the, admittance 
y^iy)) of the electromagnetic environment .Ej In terms of 
the Matsubara frequences, this relation is simply 



RKyA-^\^i\)\^i\ 

47r 



(19) 



Both S-io^ in (n5f) and the contribution to action from 
the fluctuations around classical paths are independent 
of if. Therefore, they will generate an irrelevant prefactor 
in the partition function Z. We can neglect this factor 
and obtain from (|l^ the final result 



with the total action 



(20) 



(21) 



We should point out that although So\'^ is a periodic 
function of (/?, the total action S'[(/3] is not. The physical 
origin of the non-periodic feature of S'[(^] is the contin- 
uous charge transfer between the two electrodes via the 
external circuit, which suppresses the discrete nature of 
the change of charge due to the tunneling processes. It is 
this modification of the system geometry that can smear 
out dramatically the Coulomb charging effect, as will be 
shown in the following sections. 



III. QUANTUM ADMITTANCE 

The quantum admittance i^(w) of a single tunnel junc- 
tion under a general electromagnetic environment can be 
calculated conveniently with the path integral apasciachi 
by generalizing the Kubo formula for this systenJ13^E3"E2l 



Y{lo) 



lim / dTe''^'^(T^/(r)/(0)) 



(22) 



where Tr is the time-ordering operator in imaginary time. 
The correlation function is calculated as 



(/(t)/(t')) = Z-^ J D(pe~^^-^^2e^at{T - r') 

X COs[ip{T) - ip{T')]+lT[ip,T]lT[ip,T']}, (23) 



with the current functional 



lT[^,T] = -2e dT'at{T~T')sm[^{T)~^{T')]. (24) 
Jo 

At very low temperature one can use the instanton- 
like or the renormalizatiaa .grpijp techniques to evaluate 
the partition functioi£3"EZl-E]E3, but not the correlation 
function nor the quantum conductance, which are usually 
calculated numerically.EHI However, if the temperature is 
not very low, with a semiclassical approximation we can 
evaluate analytically all path integrals to obtain a sys- 
tematic treatment of the influence of fluctuations in the 
form of a series expansion of PEc- Another advantage of 
semiclassical method is that the result is stable, Jajcause 
the fluctuation modes are of large eigenvalues .E3 Such 
relatively high temperature semiclassical results of ultra- 
small tunnel junctions turn out to be very meaningful, 
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in view of the recent measurements of Coulomb, blpckade, 
effects performed at not very low temperature.E2rLilO'L3 
In the following theoretical analysis, we will use the semi- 
classical method. 

Let us at first evaluate the partition function. The 
total action is now approximated by a Gaussian form 



dr 







l f dT f dT'a^{T-T')[^{T)-v{T')] 

^ Jo Jo 



'M2 



where 



a^ir) = at(r) + q;ox(t). 



(25) 



(26) 



The above approximation means that the contribution of 
the tunneling resistance in the action is replaced by an 
equivalent Ohmic resistance. Higher order variations of 
the action beyond the Gaussian approximataion can be 
calculated in the form of a power series of f3Ec. However, 
when calculating the current correlation function, which 
is normalized by Z, the power series correction factor in 
the denominator Z is cancelled by similar factor in the 
numerator. Hence, 5861111 given by Eq. ( psf ) turns out 
to be the approprite form for calculating the quantum ad- 
mittance as a series of PEc- We notice that an arbitrary 
path ip{T) may be expressed as a Fourier series 



(27) 



^ — — oo 



with the complex amplitudes 



ipi = ip'i+ iifi 



under the condition ip^i — (p^. Hence, the semiclassical 
action may be evaluated accordingly to give 



Sat 



1=1 



where the eigenvalues A; are 



(28) 



(29) 



Now the partition function can be directly calculated ac- 
cording to the formula 



1=1 ' 



(30) 



The current auto-correlation function ( 23 ) will be eval- 
uated in the same way. Let {I{'t)I(t'))i represent the 
first term at the right hand side of (^3|) . By performing a 
series expansion in powers of [ip{T)-ip{T')], we can express 
(J(r)J(r'))i as 



/(r)/(r'))i = 2e'at{r - t')Z~1, J D^e 



1 - 



(31) 



Making use of the Fourier series of paths '^{t), we obtain 
(/(t)/(t'))i -2e2a,(r-T') 



El „ v-^ coscjnfr — r ) 

n— 1 n— 1 



+0[{I3E,)\ 



(32) 



Substituting this expression of (/(r)/(r'))i into (p2D, we 
obtain the contribution to the admittance from the first 
term at the right hand side of (02 



Yi{u) = 2e^uj-^\ lim 



n=l 

+0[{PE, 



(33) 



The second term on the right part of ( |23| ) can be evalu- 
ated in the same way, although the algebra is more com- 
plicated. In terms of the eigenvalues of the fluctuation 
modes and the Fourier tramsform of the damping kernel, 
its contribution to the admittance is derived as 



Y2{ijj) ^ Ae^uj^'^ I lim 

iuji — >uj-{-i5 



2at{uJi) ^ atjuJi + i^n) + at{uJi - w„) - 2Q!f (tj„) 



A; 

v3/ 



n=l 

iaHuji) ^ 1 
A2 ^/3A„ 

* n— 1 



(34) 



So far our theoretical analysis is for a general elec- 
tromagnetic environment. In the rest of this paper we 
will present a detail quantitative study on dc conduc- 
tance G=lim;^^o{Iin^(^)} for a specific case that the 
electromagnetic environment contains only an Ohmic re- 
sistance. Such elec trorn agnetic environment is relevant to 
recent experimentsE2fEd with which our theoretical result 
can be compared. 
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IV. SIMPLE ELECTROMAGNETIC 
ENVIRONMENT 

We will investegate the dc conductance of a single 
tunnel junction connected in series with a resistance 
Rex- In real sample this external resistance can be var- 
ied in a controllable way. In this case we then have 
ye{—i\uJi\) = l/i?cx, and thus becomes linear in fre- 
quency 



aox(w;) 



An ■ 



(35) 



where acx—Rx /Rex is the dimensionless conductance 
contributed by the electromagnetic environment. Substi- 
tuting into ( [2S| ) this Q:cx(^i) and the at{wi) given by (||), 
we obtain the eigenvelue A; for the simple resistance en- 
vironment. If I is very large, A is dominated by the term 
ujf /2Ec, and is insensitive to the tunneling terms repre- 
sented by aex(w;) and at{uji). Hence, the corresponding 
effect turns out to be an uninteresting prefactor and con- 
tributes nothing to the conductance. For uji much less 
than the bandwidth D, the eigenvalues reduce to 



A, = 



2Er 



(36) 



where Uw—oit+acx- Inserting the eigenvalues into (33) 
after a straightforward algebra, we obtain 



Gi = lim{ImYi(w)} 



2e^ 



^J-t'^''i^^w + 1) 

2tt 



-[vI/(^„-M)-vI,(l) 



An 



(37) 



where "^{z) and "^'(z) are, respectively, the digamma and 
trigamma functions. The parameters fit and fi^ are de- 
fined as 



fXi = at(3Ec/2n^ 



t, w. 



(38) 



Now we calculate from ( pi] ) the second part of the dc 
conductance G'2=lim^_^o{Iml2('-^)}- When carrying out 
the limits at the right hand side of (34), the orders in the 
series of l3Ec become mixed. This makes the calculation 
of G2 very complicated, but yet straightforwards. The 
final result is 



Go 



+ 



Rk au 
4a? 



{aw - 

3 



RkO: 
0[{f}E,f 



-f 1) - vl/(i)]} 



*(!)] 



(39) 



Consequently, the dc conductance of a single tunneling 
junction under the resistance environment is derived as 



G = 



1 



1 - 



l3Ec ae 



Rex + Rt 
^(/i^ + 1) - ^(1) 



«''(/i. 



(40) 



In the high temperature limit, the tunneling system ex- 
hibits the conventional Ohmic behavior Gohm — 1/ {Rt + 
Rex) as expected. 



V. RESULTS AND DISCUSSION 

To demonstrate the influence of electromagnetic envi- 
ronment on Coulomb charging effects, let us introduce 
the dimensionless conductance correction 



Ag = l- G{R, 



Ri 



_ PEc ag: 



*(Ai^ + 1) - *(i) 



+ i''(^„ + 1) 



(41) 



which can be considered as a measure of the effective 
Coulomb charging energy: smaller corresponds to 
weaker effective. Coulomb charging energy. In recent 
experimentaH30 Pekola and his co-workers have man- 
ufactured a single tunnel junction with at=i?^/i?t=0.5 
and measured I/A5 as a function of normalized temper- 
ature ksT/Ec for different values of acx—RK/Rcx- Rex 



contains two contributions Rc.x=Ro 



Re 



where 



-Rex, in is the intrinsic value of the circuit, which is of the 
order of the free space impedance \/ Ho/eo — 377 Q and 
J^cx.con is experimentally controlable. At high tempera- 
ture 1 / Ag is found to be linear in ksT/ Ec- Extrapolating 
this high temperature straight line to zero temperature, 
it was found that the zero temperature offset {1/ /S.g)Q is 
positive and decreases with diminishing aox- Our theory 
has reproduced these observed features as shown in Fig. 1 
for at—RK/Rt—0-^- For high temperature, 1/Ag can be 
readily derived as 



1 Saw ksT 
Ag Qfcx Ec 



(0.6 + 0.2a„)a« 



(42) 



Thus, the offset is simply (l/A5)o = (0.6 -|- 
0.2aw)c(w/ctcx- For the case that the tunnel conductance 
at is much smaller than the electromagnetic environ- 
mental conductance aox, and in the vicinity of fiw'til^ 
the offset can be well approximated by the very sim- 
ple form (1/A(7)o=0.6 -I- 0.2Rk / Rex- Hence, if we take 

-Rex, in — 377 we obtain (l/A(/)o — 14 for -Rex, con =0 and 

(l/Ag)o — 2 for i?exp™rrj=3 kil, which agree exactly with 
the observed valuesxStJ An effective-capacitance model 
has been proposecl23'Ell to explain these measured offset 
values. The so-calculated {l/Ag)Q is smaller than the 
experimental value by a factor 5 for -Rox,con=0, and fails 
to explain the case of -Rex,con=3 kQ. 
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For a fixed value of PEc, Aij is a function of the junc- 
tion conductance at and the electromagnetic environ- 
mental conductanGe-|acx- Using a standard formalism, 
it has been proverO that for high temperature, in the 
plane of at and aox, the ratio Ag/ f3Ec has an absolute 
maximum value 1/3. However, the functional depen- 
dence of Ag/PEc on at and aox has never been obtained. 
Our theory allows us to calculate such important results 
which are shown in Fig. 2 for /3£'c=0.2 and in Fig. 3 for 
PEc=0.02. In each figure the ratio n^Ag/PEc is plotted 
as a function of aex for at=0.1 (solid curve), 1.0 (dot- 
dashed curve), 10.0 (dotted curve), and 100.0 (dashed 
curve). It is clear that in the at-acx plane, the abso- 
lute maximum value of Ag / PEc lies in the region where 
Rt^Rcx> Rk- a careful numerical search reveals that 
the value of this absolute maximum is just 1/3, corre- 
sponding to the charging effects due to the geometric 
capacitance of the ultrasmall tunnel junction. 

We notices that in both Fig. 2 and Fig. 3, for a fixed 
value of acx, the value of ir-^Ag/PEc decreases monotoni- 
cally with at ■ This is understandable because an increase 
of the junction conductance (or the tunnel strength) leads 
to a reduction of the effective Coulomb charging energy. 
On the other hand, for a fixed value of the junction con- 
ductance at, n'^Ag/PEc as a function of the external con- 
ductance acx exhibits a peak. Such interesting feature 
has been obsepfed very recently, but has not appeared in 
literature yet.Eil The physics of this peak structure can be 
explained as follow. Once an electron tunnels through the 
junction from the left side to the right side, it can return 
to the left side via the external circuit. Therefore, the 
larger is the external conductance aox, the faster can the 
tunnel junction relax to its initial charge configuration. 
Consequently, the effective Coulomb charging energy is 
reduced and so is Ag. This conclusion then agrees with 
the one reached in Ref. |^ and Ref. ||. It is worthwhile to 
mention that in the limit of very high tunneling rate or 
very large aox such that fi^^^l, (^) is simplified to 



Ag. 



2ao 



ln{ 



OwPEc 

27r2 



} + l-vE.(l) 



(43) 



which is not analytical in at. Therefore, the method fo 
calculating G in Ref. ^ based on the Fermi golden rule 
is no longer valid here. 

The picture is entirely different at the other end of very 
small external conductance aox, or very large external 
resistance Rcx/Rk- In this case the external returning 
path from the right side of the junction to the left side is 
heavily blocked. Hence, once an electron tunnels through 
the junction, its probability to tunnel back is much larger 
than that to travel through the external circuit. Eventu- 
ally, an electron spends most of the time tunneling back 
and forth, and so higher order tunneling terms must be 
taken into account, even though each tunnel process is 
still incoherently sequential. Such tunneling processes 
cannot be treated with the Fermi golden rule neither. In 
fact, the system is in a dynamical state with an equivalent 



reduced effective Coulomb charging energy. Each curve 
in Fig. 2 and Fig. 3 then drops down to zero as aex ap- 
proaches zero, at which the external circuit is completely 
blocked. This point of view is similar to that appeares 
in the investigatipp-joLthe statistical behavior of a single 
electron box.EElEJ Nevertheless, in the present prob- 
lem the renormalization of the capacitance is dynamical 
rather than static. To our knowledge, such dynamical 
renormalization of the charging energy in the regime of 
very small external conductance aox is discovered for the 
first time. 



VI. REMARKS 

We have used a non-perturbative approach to derive 
the quantum admittance of a single tunnel junction in the 
presence of an electromagnetic environment. Besides the 
case with an Ohmic environment which has been inves- 
tigated in the present work, our theoretical formulation 
is very general and can be conveniently used to study 
the effects of non- Ohmic environments. If the electro- 
magnetic environment contains a frenij.PT].ry-dependent 
part described by a transmission lineO'Lfcjtj, the quan- 
tum admittance of the system as a function of frequency, 
which is of great interesting, can be analyzed with our 
theory. Such theoretical results can iift-compared with 
experiments in a very direct manner.EjC3 With a slight 
modification of the action, we can also investigate the 
behavior of an array of tunnel junctions. If the junction 
array has a pure Ohmic environment, the dc conductance 
of the system can be calculated by evaluating the piean 
value of the current through the external resistance.EZI Fi- 
nally, we should mention that very recently a dissipation- 
driven superconductor-insulator transition in a Joseph- 
son junction array has been observed by changing the 
external impedance continuously.ES Using the theoretical 
approach developed in the present work, this problem 
can be investigated theoretically at a quantitative level. 
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FIG. 1. 1/Ag as a function of the normalized temperature 
ksT/Ec for the dimensionless junction conductance at=0.5, 
with various values of the dimensionless external conductance 
aex=10 (solid curve), 50 (dotted curve), 100 (dashed curve), 
and 150 (dot-dashed curve). The cases Qox=10 and 100 cor- 
respond closely to the two measured samples in Ref. 20 and 
Ref. 21, respectively. 



FIG. 2. Ag in units of PEc/n^ as a function of the di- 
mensionless external conductance ttex for f3Ec=0.2, with var- 
ious values of the dimensionless junction conductance at=0.1 
(solid curve), 1.0 (dot-dashed curve), 10.0 (dotted curve) and 
100.0 (dashed curve). 



FIG. 3. The same as Fig. 2 but for PEc=0.02. 
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